AD-A107  669 

UNCLASSIFIED 


GEORGE  WASHINGTON  UNIV  WASHINGTON  DC  PROGRAM  IN  LOGISTICS  F/6  13/13 
A  CLOSED  FORM  LOCAL  SOLUTION  OF  A  NONLINEAR  STRUCTURAL  DESIGN  P—ETC(U) 
AUG  01  A  V  FIACCOt  A  GHAEMI  N00014— 75— r  — 0729 

SERlAL-T-449  n. 


STUDENTS  FACULTY  STUDY  R 
ESEARCH  DEVELOPMENT  FU7 
URE  CAREER  CREATIVITY  CC 
MMUNITY  LEADERSHIP  TECh 
NOLOGY  FRONTIE 
ENGINEERING  APfife^lBIE 
GEORGE  WASHIh 


DTIC 

ELECTE| 
NOV  20  1981 


SCHOOL  OF  ENGINEERING 
AND  APPLIED  SCIENCE 


C  ^  •  t  J72 


mu  document  has  kcn  Armovco  rot  wtuc  ieuak  and  sau;  its  outkikition  is  unumiteo 


A  CLOSED  FORM  LOCAL  SOLUTION  OF  A 
NONLINEAR  STRUCTURAL  DESIGN  PROBLEM  IN 
TERMS  OF  TEE  DESIGN  PARAMETERS 


by 

Anthony  V.  Fiacco 
Abolfazl  Ghaemi 


Serial  T-449 
10  August  1981 


The  George  Washington  University 
School  of  Engineering  and  Applied  Science 
Institute  for  Management  Science  and  Engineering 


Program  iri  Logistics 


Accession  For 

NT  I S  GRA&I 
DTIC  TAB 
Unannounced 


Justification 


By - - — 

^Distribution/  _ 
Availability  Codes 
jAvaii  and/or 
Diet  I  Special 


ft 


Contract  N00014-75-C-0729 
Project  NR  347  020 
Office  of  Naval  Research 


This  document  has  been  approved  for  public 
sale  and  release;  its  distribution  is  unlimited* 


DTIC 

SELECTED 
NOV  20  19811  ■ 


□  □ 


SECURITY  CLASSIFICATION  OF  THIS  FAGC  f*h*n  Dete  Entered) 


REPORT  DOCUMENTATION  PAGE 


RtPQRT  NUMBER 


^  T-449 _ 


4.  title 

A  CLOSED  FORM  LOCAL  SOLUTION  OF  A 
'  NONLINEAR  STRUCTURAL  DESIGN  PROBLEM  IN 
TERMS  OF  THE  DESIGN  PARAMETERS  - 


7.  author^; 

ANTHONY  V./FIACCO 
ABOLFAZL  GHAEMI 


»  PERFORMING  ORGANIZATION  NAME  AND  ADORES* 

THE  GEORGE  WASHINGTON  UNIVERSITY 
PROGRAM  IN  LOGISTICS 
WASHINGTON,  DC  20052 


11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

OFFICE  OF  NAVAL  RESEARCH 
CODE  434 

ARLINGTON,  VA  22217 


MONITORING  AGENCY  NAME  *  ADDRESS^//  different  from  Controlling  Oftice J 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


1  RECIPIENT'S  CATALOG  NUMBER 


»  TYPE  OF  REPORT  *  PERIOD  COVERED 


kTHIDUHSIJIH 


«.  .crpojwuuj  one.  r«rort  nuukk 

T-449 


•  contract  or  grant  num.ero; 


N00014-75-C-0729 


IS.  REPORT  DATE 

10  AUGUST  1981 


1).  NUMBER  OP  PAGES 

22 


Controlling  Office)  IS-  SECURITY  CLASS,  (of  this  report) 

UNCLASSIFIED 

DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 


IS  DISTRIBUTION  STATEMENT  (of  thla  Haport) 


APPROVED  FOR  PUBLIC  SALE  AND  RELEASE;  DISTRIBUTION  UNLIMITED. 


17.  DISTRIBUTION  STATEMENT  (of  tho  ebatrect  entered  In  Block  20,  It  different  from  Heperl) 


IS-  KEY  WOROS  ( Continue  on  reveree  oido  If  neceeeery  end  Identify  by  block  number) 

NONLINEAR  PROGRAMMING 
STRUCTURAL  DESIGN 
PARAMETERS 

CLOSED  FORM  SOLUTION 


20.  ABSTRACT  (Continue  on  reveree  aide  It  neceeeery  end  Identity  by  block  number) 

\  _  _ _ _ _ _ _ 

Using  the  results  of  our  previous  analysis  of  the  corrugated  bulkhead 
model,  a  closed  form  of  the  optimal  solution  of  this  model  is  easily  de¬ 
rived  as  a  function  of  the  many  design  parameters.  The  analytic  solution 
is  valid  over  large  changes  in  the  parameters.  This  demonstrates  the 
possibility  of  making  very  accurate  (in  this  instance,  exact)  local  estimates 
of  a  markedly  nonlinear  parametric  solution,  once  one  solution  is  found, 


DD  1473 


EDITION  or  1  NOV  «i  IS  OBSOLETE 
S/N  0103*014-4601 


(continued  ) 


_ UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Dote  Stored) 


IHITV  CLASSIFICATION  or  THIS  PA OCfPh«w  Dtm  fhMMtf; 


ABSTRACT  (cont'd) 


/  even  though  the  original  problem  statement  might  appear  to  make  this 
prohibitive.  It  is  shown  that  the  local  parametric  solution  obtained 
can  be  extended  to  include  large  changes  in  the  parameters ,  particularly 
when  the  parametric  perturbations  are  restricted  and  highly  structured, 
as  is  typical  in  practical  applications.  Aside  from  the  practical  value 
of  such  a  result,  it  is  reported  because  this  problem  is  well  publicized 
and  has  been  utilized  for  some  time  as  a  test  problem.  Knowledge  of  the 
precise  solution  for  a  range  of  parameter  values  provides  Ideal  information 
for  further  studies.  , 

r 


•ceufttrv  ciMtiricATioN  or  this  pa  sap***  d«h 


if 


THE  GEORGE  WASHINGTON  UNIVERSITY 
School  of  Engineering  and  Applied  Science 
Institute  for  Management  Science  and  Engineering 

Program  in  Logistics 


Abstract 

of 

Serial  T-449 
10  August  1981 


A  CLOSED  FORM  LOCAL  SOLUTION  OF  A 
NONLINEAR  STRUCTURAL  DESIGN  PROBLEM  IN 
TERMS  OF  THE  DESIGN  PARAMETERS 


by 

Anthony  V.  Fiacco 
Abolfazl  Ghaemi 


Using  the  results  of  our  previous  analysis  of  the  corrugated 
bulkhead  model,  a  closed  form  of  the  optimal  solution  of  this  model  is 
easily  derived  as  a  function  of  the  many  design  parameters*  The  analytic 
solution  is  valid  over  large  changes  in  the  parameters.  This  demonstrates 
the  possibility  of  making  very  accurate  (in  this  Instance,  exact)  local 
estimates  of  a  markedly  nonlinear  parametric  solution,  once  one  solution 
is  found,  even  though  the  original  problem  statement  might  appear  to  make 
this  prohibitive.  It  is  shown  that  the  local  parametric  solution  obtained 
can  be  extended  to  include  large  changes  in  the  parameters,  particularly 
when  the  parametric  perturbations  are  restricted  and  highly  structured, 
as  is  typical  in  practical  applications.  Aside  from  the  practical 
value  of  such  a  result,  it  is  reported  because  this  problem  is  well 
publicized  and  has  been  utilized  for  some  time  as  a  test  problem. 

Knowledge  of  the  precise  solution  for  a  range  of  parameter  values  pro* 
vides  ideal  information  for  further  studies. 


Research  Supported  by 

Contract  N00014-75-C-0729 
Project  NR  347  020 
Office  of  Naval  Research 


i 


TABLE  OF  CONTENTS 


Abstract  i 

1.  THE  PROBLEM  1 

2.  DERIVATION  OF  CLOSED  FORM  SOLUTION  11 

3.  NUMERICAL  CORROBORATION  14 

4.  DEVELOPMENT  OF  A  COMPLETE  PARAMETER  SOLUTION  16 

REFERENCES  19 


t 

t 


il 


E 


i 


THE  GEORGE  WASHINGTON  UNIVERSITY 


School  of  Engineering  and  Applied  Science 
Institute  for  Management  Science  and  Engineering 


A  CLOSED  FORM  LOCAL  SOLUTION  OF  A 
NONLINEAR  STRUCTURAL  DESIGN  PROBLEM  IN 
TERMS  OF  THE  DESIGN  PARAMETERS 


by 

Anthony  V.  Fiacco 
Abolfazl  Ghaemi 


1.  THE  PROBLEM 


The  purpose  of  this  study  is  to  further  characterize  the  optimal 
solution  behavior  of  a  nonlinear  programming  model  for  the  optimal  sizing 
of  a  vertically  corrugated  transverse  bulkhead  of  an  oil  tanker ,  presented 
in  [l,Ch.  6]  and  subsequently  analyzed  by  the  authors  [2].  The  objective 
of  the  model  is  to  determine  the  dimensions  of  the  bulkhead  such  that 
certain  engineering  constraints  are  met,  with  minimum  possible  bulkhead 
weight. 

Figures  1  and  2  depict  the  configuration  of  the  bulkhead  and 
design  variables,  respectively.  Figure  3,  together  with  Table  1,  defines 
the  model  parameters  and  their  values. 
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q-Yhlb«  (load) 


Figure  3.  Specification  of  some  of  the  design  parameters 
and  lndlcatidn  of  load  levels,  side  view. 


TABLE  1 


PROBLEM  DESIGN  PARAMETERS 
(MODEL  INPUT  DATA) 


Parameter 


Meaning 


Weight  per  unit  volume  of  the 
material 


Width  of  the  panel 

Length  of  the  top  panel 

Length  of  the  middle  panel 

Length  of  the  bottom  panel 

Distance  between  free  liquid 
level  and  top  of  structure 

Distance  between  free  liquid 
level  and  middle  of  top  panel 

Distance  between  free  liquid 
level  and  middle  of  middle 
panel 

Distance  between  free  liquid 
level  and  middle  of  bottom 
panel 

Distance  between  free  liquid 
level  and  base  of  top  panel 

Distance  between  free  liquid 
level  and  base  of  middle  panel 

Distance  between  free  liquid 
level  and  base  of  bottom  panel 

Minimum  allowable  thickness  of 
top  panel 


476  cm 
495  cm 


385  cm 


315  cm 


250  cm 


498  cm 


938  cm 


1288  cm 


745  cm 


1130  cm 


1445  cm 


1.05  cm 


Minimum  allowable  thickness  of 
middle  panel 


1.05  cm 
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Table  1  —  Continued 


15 


Parameter 


min 

Cb 


Meaning 


Minimum  allowable  thickness  of 
bottom  panel 


Value 


1.05  cm 


16 


e 


Effectiveness  of  the  flange 
(dimensionless) 


0.8 


17 


A  derived  coefficient  (func¬ 
tion  of  maximum  allowable 
bending  stress) 


6.94  x  10 


-8 


-1 

cm 


18 


Corrosion  coefficient 


0.15  cm 


The  resulting  nonlinear  programming  problem  obtained  in  [1]  takes 
the  following  form. 

minx  f  =  rB(xx  +  x2)(£tx4  +  Axs  +  Z^)  ^  +  (x2  -  x2)172]"1  (weight) 


subject  to 

(i)  Geometrical  constraint: 
gl  5  x2  '  x3  =  °  ; 


-  6  - 
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(il)  Bending  stress  constraints: 

2  22  1/2 

g2  H  X2X3X4  +  3e  X1X3X4  “  ^klht*t  ^X1  +  ^x2  “  x3^  1  i  0  ; 


g3 


x_x,xc  +  3e  x,x,xc  -  6k 


.  02  r  ,  ,  2  2.1/2  s  . 
lhm  m  [xl  +  X2  "  X3J  -  0  5 


2  22  1/2 

g4  =  *2*3*6  +  3e  X1X3X6  ~  6kihb^b  ^X1  +  ^x2  “  *3^  ]  —  0  ; 


(iii)  Moment  of  inertia  constraints 
2 


g5  =  X2X3X4 


S6  “  X2X3X5 


g?  =  x2x3x6 


.  ,  2  ?2. 4/3,  .  ,  2  2. 1/2.4/3  *  A 

+  3e  XjX^x^  -  26.4(k1ht*-t)  ^  +  (x2  -  x3)  ]  >0 

,  ,  2  ,,  ...  .  «2,4/3r  .  ,  2  2* 1/2.4/3  ^  A 

+  3e  x.x,xc  -  26.4(k,h  * _)  [x,  +  (x_  -  x,)  ]  >.  0 

.  ,  2  ......  »2.4/3r  .  ,  2  2. 1/2.4/3  .  A 

+  3e  x^x^  -  26.4(k1hb*.b)  [x^  +  (x2  -  ,:3>  ]  >  0 


(iv)  Thickness  requirement  constraints: 


.min  .  A 

g8  =  x4  “  Ct  =  0  ; 


1/2, 


gg  =  x4  -  [.39  •  1.05(.01h1J;)i'‘-(.01x1)  +  k2]  >  0  ; 
glQ  H  x4  -  [.39  •  1.05(.01hlt)1/2(.01x2)  +  k2]  >  0  ; 

811  5  X5  '  i  0  5 

1/2 

g12  =  x5  -  [.39  •  1.05(.01hlm)x/  (.01Xl)  +  k2]  >  0  ; 

1/2 

gl3  =  x5  -  [.39  •  1.05(.01hlm)1/z(.01x2)  +  k2]  >  0  ; 

min 


g  =  x6  ~  >  0  ; 


,1/2, 


g15  5  x6  "  1,39  '  l-05C.0Uilb)- ^‘(.Olx^  +  k2]  >  0  ; 

1/2 

g16  =  x6  "  1,39  ‘  l-05(.01hlb)x/^(.01x2)  +  k2]  >  0  ; 


(v)  Natural  constraints: 


g17  *  g18:  xi  ^  0  1  "  1*3 


-  7  - 
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Table  2  gives  the  optimal  solution  and  Table  3  shows  the  constraint 
values  and  Lagrange  multipliers  corresponding  to  this  solution.  These 
results  were  obtained  in  our  previous  study  [1]. 


TABLE  2 

OPTIMAL  SOLUTIONS  OF  THE  MODEL 
WITH  PARAMETER  VALUES  GIVEN  IN  TABLE  1 


X1 

* 

bl 

= 

57.82 

cm 

X2 

b2 

= 

57.82 

M 

X3 

- 

d 

m 

35.69 

tt 

X, 

3 

t 

3 

1.05 

n 

4 

t 

xc 

t 

B 

1.05 

H 

5 

m 

x6 

= 

cb 

m 

1.05 

ii 

f 

m 

w 

s 

5.25 

tons 
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TABLE  3 


LAGRANGE  MULTIPLIERS  AND  CONSTRAINT  VALUES  AT  THE  OPTIMAL  SOLUTION 
(excluding  nonnegativity  constraints) 


Constraint 

8i 

Value  of 
Constraint 

Value  of  Lagrange 
Multiplier 

81 

, 

22.13 

.451878-*  x  10~2 

*2 

2117.22 

.4723479  x  10"4 

8^ 

1383.41 

.7218857  x  10"4 

*4 

1868.92 

.5351946  x  10”4 

*5 

41993.15 

.2381633  x  10"5 

.043 

.2300013  x  101 

*7 

27945.57 

.3580375  x  10“5 

.48  *  10"5 

.2078305  x  107 

*9 

.25 

.3940568 

g10 

.25 

.3940570 

+811 

.96  x  10~5 

.1040634  x  107 

*12 

.10 

.9604345 

*13 

.10 

.9604360 

+814 

.49  x  iff4 

.2029215  x  106 

.31  x  l<f4 

.3221582  x  106 

+816 

.13  x  i<f 4 

.7974950  x  106 

tBlnding  constraint* 
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As  indicated  in  Table  3,  constraints  »  8g  »  gj^  and 

gw  -  g,,  are  binding.  It  was  shown  in  (2]  that  the  corresponding 
14  lo 

gradients  of  these  functions  are  linearly  independent.  The  optimal 
solution  of  the  model,  with  the  given  data  base  (Table  1),  thus  corre¬ 
sponds  to  the  locally  unique  solution  of  the  following  nondegenerate 
system  of  nonlinear  equations: 


2  A  ,  2  „  p  2v 4/3r  .  ,  2  2. 1/2,473  A 

g6  *  X2X3X5  +  3e  X1X3X5  *  26,4^kl  m  ^X1  +  ^*2  “  *3'  J  *  0 


8 

*  x,  - 

tmin  =  0 

4 

t 

11 

*  X-  — 

tmin  =  0 

5 

m 

14 

=  X6‘ 

t”in  =  0 

b 

15 

=  X6  “ 

1/2 

£(.39)(1.05)(.01hlbr'  (.Olx^  +  k2]  *  0 

16 

*  X6  " 

[(.39)(1.05)(.01hlb)1/2(.01x2)  +  k2]  -  0 

In  [2]  we  determined  that,  with  parameter  perturbations  as  large 
as  two  percent  of  the  parameter  base  values,  the  set  of  bending  constraints 
remains  unchanged.  This  means  that,  for  relatively  large  changes  of  data 
in  a  neighborhood  of  the  data  base  given  in  Table  1,  the  optimal  para¬ 
metric  solution  of  the  nonlinear  program  is  unique  and  must  satisfy  the 
system  of  equation  (1). 

In  the  next  section,  ve  obtain  a  closed  form  solution  for  this 
system  of  nonlinear  equations. 


-10  - 
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2.  DERIVATION  OF  CLOSED  FORM  SOLUTION 


For  notational  simplicity,  we  denote  the  coefficents  of  Equations 


(1)  as  follows: 


C1  =  3e 


c.  -  26.4(k,h  £2)4/3 

2  1  ra  m 


c3  “  k2 


min 

C5  *  Sn 


_min 

c6  '  S 


c,  -  (.39)(1.05)(.01hlb)  '  (.01) 


Thus,  the  system  of  equation  (1)  becomes 


2  2  r  .  ,  2  2,1/2,473 

g6  :  +  c1x1x3x3  -  c2^xl  +  ^x2  ~  *3)  *  *  0 


g8  :  x4  ‘  c4  *  0 


gll  :  x5  "  c5  "  0 


g14  :  x6  ‘  c6  *  0 


g15  :  x6  -  C7X1  "  C3  "  0 


g16  :  x6  ~  C7X2  “  c3  “  0 
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From 


88  ’  gll  3nd  814  W<i  obtain 


x,  =  c.  ,  x,  ■  c,  and  x,  -  c, 

4  4  5  5  6  6 


(3) 


From 


g15  and  g16  and  (3)  , 


xx  «  x2  =  (c6  -  c3)/c7  =  Cg 


(4) 


Substituting  x^  and  Irom  (4)  and  x,.  from  (3)  in  g^  ,  we  obtain 

/  ^  ,2  r  ,  ✓  2  2,l/2,4/3  . 

(C5C8  +  C1C8C5)  X3  "  C2  C8  +  (c8  _  X3}  1  *  0 


or 


,  ^  j.  ,  /  ,3/4  3/2  r  .  ,  2  2,1/2, 

(c5cg(l  +  c1)/c2^  x3  *  lc8  +  (cg  -  X3)  ]  , 


/  csc«(1  +  CP\3^4 
Letting  c,  -  (-5-2- - 2)  , 


(5) 


then  subtracting  Cg  from  both  sides  of  the  equation  and  squaring. 


yields 


or 


2  3  2  0  3/2 

c9x3  +  x3  -  2c8c9x3  -  0 


3/2f  2  3/2  .  1/2  ,  , 

x3  [c9  x3  +  x3  "  2c8C9 


Since  we  must  have  >  0  near  the  given  base  value  solution  and 
since  we  may  assume  >  c5  >  c2  >  ^  atw*  c6  >  c3  *  then  >  0 


and  Cg  +  0  ,  so 


8 


(TxT?  +  ~2  7x7  -  — 

3  Cy  3  9 


reduced  cubic  equation  in  7x7  . 


(6) 
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According  to  a  classical  result  discovered  by  Ferro  and  first 
published  by  Cardan  [3],  a  reduced  cubic  equation  of  the  form 

+  ay  +  b  *  0  (7) 

with  (y)^  +  (y) ^  >  0  has  only  one  real  root. 

Since  in  (6)  ,  a  *  — y—  is  always  positive,  (yj^  +  (y)~*  is 

C9 

always  positive  for  our  problem  for  all  changes  of  these  coefficients 
that  result  from  any  (feasible)  model  parameter  perturbations. 

The  unique  real  root  of  (7)  is  given  by  Cardan's  formula  as 
follows: 


thus,  the  solution  for  x  from  (6)  is 


We  have  thus  obtained  a  closed  form  parametric  solution  for  the 
optimal  problem  variables  near  the  given  initial  parameter  values. 
Substituting  -  x^  from  (3),  (4)  and  (8)  in  the  objective  function 

f  of  the  model  well  of  course,  yield  the  closed  form  formula  for  the 
optimal  value  function. 
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3 .  NUMERICAL  CORROBORATION 


The  closed  formulas  we  have  derived  provide  an  opportunity  to 
check  the  computer  solution.  In  this  section,  we  use  the  data  base 
given  in  Table  1  and  the  solution  formulas  we  have  derived  to  calculate 
the  optimal  solution  components  and  objective  function  value,  and  compare 
these  with  those  obtained  by  the  computer  solution  and  listed  in  Table  2. 

Using  the  data  in  Table  1,  we  obtain  the  following  values  for 


and 


and 


C1  *  2,4  *  c2  *  542. 311667  ,  c3  -  .15  •  c4  “  c5  *  c6  "  1-05 
c  ?  *  .015566  . 


From 

<«>  , 

c  -  c 

c„  -  — - i  »  57.818322  , 

8  c7 

from 

(5)  , 

/c  c  (1+c  )\3/4 

C9-l  c2  )  "  *484578 

Substitution  of  c^  and  c^  in  (8)  yields. 

x3  -  (5.973894)2  -  35.687409  . 


Thus,  the  analytically  calculated  solution  vector  is,  to  two  decimal 
points  of  accuracy. 


x  -  (57.82  ,  57.82  ,  35.69  ,  1.05  ,  1.05  ,  1.05)  , 

which  corresponds  precisely  to  the  computer  solution  given  in  Table  2. 
Substantiation  of  these  values  of  x  in  the  objective  function  f  and 
rounding  to  two  decimal  points  yields  a  value  of  5.25  ,  again  in  precise 
correspondence  with  the  computer  solution. 
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As  an  example  of  a  verification  of  previously  derived  sensitivity 
results,  we  calculate  the  sensitivity  of  the  optimal  value  function  to 

the  parameter  t^n  ,  the  right-hand  side  of  the  constraint  gg  which 

we  denoted  by  • 

Denoting  x^  by  and  substituting  the  closed  form  solutions 

for  x^  i=l,6  in  the  objective  function  f  ,  we  obtain  the  formula 
for  the  optimal  value  function. 


* 

f 


(•) 


2rBc8(t'tC4  + 


;5  +  W 


2  2 
(c8  ”  C10) 


1/2 


+  cr 


(9) 


Thus, 


9c4 


2rBc8At 


c2  )1/2 
10' 


Substituting  in  (9)  the  numerical  values 
in  Table  1  and  the  calculated  values  of 
above  formulas,  we  obtain 


of  T  ,  B  and  given 

Cg  and  c^q  given  in  the 


3fV) 

9c4 


2070337.126  . 


The  computer  solution  obtained  in  [1]  was  2078305  .  This  corresponds 
with  the  calculated  Lagrange  multipliers  value  associated  with  constraint 
gg  ,  which  is  given  in  Table  3.  These  calculations  are  in  close  agree¬ 
ment,  in  terms  of  relative  error. 
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4.  DEVELOPMENT  OF  A  COMPLETE  PARAMETRIC  SOLUTION 


A  careful  study  of  the  local  solution  that  has  been  obtained  and 
the  constraint  structure  indicates  that  a  great  deal  of  progress  can  be 
made  towards  obtaining  a  global  parametric  solution,  particularly  if 
the  parameter  perturbations  are  restricted  to  those  that  typically 
arise  in  practice.  For  example,  perturbations  of  considerable  interest 
are  based  on  the  assumption  that  some  or  all  of  the  parameters  change 
"proportionately."  e.g.,  if  a  given  parameter  vector  is  denoted  by 
p  ,  then  a  perturbed  value  of  p  would  be  given  by  p  +  ap  -  (1  +  a) p  * 

3p  ,  where  a  is  the  proportionality  factor  and  6  =  1  +  a  ,  a  scalar. 
With  this  restriction,  all  the  parameters  that  change  are  functions  of 
the  single  scalar  valued  parameter  3  and  the  perturbation  analysis 
is  enormously  simplified.  In  the  present  example,  given  this  type  of 
perturbation,  we  could  readily  track  the  course  of  the  closed  form 
parametric  solution  that  we  have  obtained  as  a  function  of  the  single 
parameter  3  ,  as  3  increases  or  decreases.  As  long  as  the  given 
binding  constraints  remain  binding,  the  given  system  of  equations  hold 
and  the  given  parametric  solution  would  apply.  When  a  new  constraint 
becomes  binding  at  a  new  value  of  3  ,  we  would  have  to  introduce  a  new 
equation  and  when  a  binding  constraint  becomes  nonbinding,  the  correspon¬ 
ding  equation  would  be  removed  from  the  binding  constraint  equations. 
Given  the  relative  simplicity  of  our  closed  form  solution  and  the 
exploitable  structure  of  the  constraints,  as  we  shall  indicate,  a 
complete  closed  form  parametric  solution  could  probably  be  developed 
over  the  domain  of  parameter  values  of  interest. 

At  a  general  level,  note  that  constraints  and  gg- 

are  simple  linear  constraints.  Each  of  these  constraints  that  are 
binding  will  either  determine  the  value  of  a  variable  or  determine  one 
variable  in  terms  of  another,  thus,  each  removing  one  degree  of  freedom 
in  the  solution.  If  any  of  the  constraints  g^  -  are  binding,  they 

can  each  be  solved  for  any  variable  in  terms  of  the  others  involved  in 
the  particular  constraint.  Finally,  if  either  of  g^  -  are  binding 
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then,  x.  ,  x,  ,  xc  ,  or  x,  ,  can  be  expressed  in  terms  of  the  other 

1  H  D  O 

variables,  and  if  x^  =  x2  >  then  the  remaining  variables,  x^  and  , 

as  we  have  shown,  can  be  solved  in  terms  of  the  other  variables.  Thus, 
the  preponderance  of  the  various  possible  combinations  of  equations  that 
can  arise  from  the  constraints  that  are  binding  at  a  solution  corresponding 
to  a  given  parameter  value  can  be  resolved  in  closed  form. 


We  conclude  with  an  extension  of  the  local  closed  solution  that 
we  have  obtained  to  illustrate  our  observations.  Consider  changing  the 
parameter  t®in  from  its  base  value  1.05  cm  ,  leaving  the  other  para¬ 
meters  fixed  as  given  in  Table  1.  Denote  the  parameter  vector  by  £  , 
so  that  denotes  the  ith  parameter  value,  and  the  closed  form 

solution  vector  that  we  have  obtained  by  x(e)  .  If  =  t”*” 

increases  from  its  base  value  to  any  larger  value  t®*n  ,  with  the  other 

* 

fixed  at  their  initial  values  ei  ,  then  the  optimal  parametric 

solution  vector  is  given  by  x^(£)  38  t™*n  9  w^t^1  the  other  components 

min 

x  (c)  of  the  solution  unchanged.  On  the  other  hand,  if  t^  decreases 
to  ,  with  the  other  parameter  values  fixed  at  their  given  initial 

values,  then  we  find  that  the  optimal  parametric  solution  is  given  by 


*A(t> 


,  with  no  change  in  all  the  other  components,  providing 


tmin  _  tmin  sroall#  However,  the  nonbinding  constraints  82  *  g5  ’  g9 
and  g  ^  evaluated  at  this  solution  decrease  so  the  solution  remains 
optimal  providing 


>  max  (m2  ,  m3  ,  ,  m1Q) 


where  is  the  quantity  involved  in  the  ith  constraint  when  the 

ith  constraint  is  equivalently  expressed  in  the  form  >  m^[x(£),e]  , 

1-2,  5,  9,  10,  and  the  are  evaluated  at  the  given  solution  x(e) 
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for  the  given  parameter  data  base  £  .  Hence,  the  given  closed  form 
parametric  solution  is  valid  for  t^*n  >  ,  where 

t™in  *  max  (m2  ,  m5  ,  mg  ,  m1Q)  (10) 

Clearly,  g6  ,  gQ  ,  gn  ,  g14  ,  g15  and  g16  will  remain  binding,  while 

g2  >  g5  *  or  g1Q  will  become  binding  when  =  t™in  . 

We  do  not  pursue  other  possible  extensions  of  the  parametric  solu¬ 
tion.  A  moral  of  the  foregoing  analysis  is  that  closed  form  parametric 

solutions  are  not  the  exclusive  province  of  "toy"  problems.  Closed  form 

solutions  provide  precise  and  complete  information,  a  goal  that  may  be 
perfunctorily  relinquished  for  "practical  problems."  Even  if  a  complete 
parametric  solution  cannot  be  obtained,  various  components  may  be  possible 
to  derive  in  closed  form  over  parameter  perturbations  of  interest,  thus 
greatly  simplifying  further  analysis.  In  this  problem,  for  example,  even 
if  x^(£)  could  not  have  been  obtained,  the  other  components  of  the  local 
parametric  solution  are  immediately  available  in  closed  form. 

Obviously,  the  existence  of  a  complete  closed  form  parametric 
solution  is  a  fortunate  happenstance.  Components  of  nonlinear  program¬ 
ming  parametric  solutions  will  inevitably  have  to  be  estimated,  perhaps 
most  aptly  by  predicator-corrector  techniques  and  their  contemporary 
refinements,  e.g.,  continuation  methods  [4],  [5].  Nonetheless,  experience 
suggests  that  some  closed  form  functional  relationships  characterizing 
a  solution  will  inevitably  be  available  and  are  worth  exploiting  once  a 
solution  has  been  determined. 
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